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ABSTRACT 

We  present  a  simple  model  for  the  structure  constant  of  atmospheric  refractive  index  fluc¬ 
tuations  (Cn2)  and  inner  scale  ( i0 )  valid  in  the  lowest  few  hundred  meters  above  land.  The 
inputs  to  the  model  are  standard  meteorological  parameters  plus  solar  irradiance  and 
parameters  characteristic  of  the  local  terrain.  Results  are  compared  with  measurements  of 
Cn  and  60  made  with  a  displaced-beam  scintillometer.  The  results  show  that  the  model  can 
predict  atmospheric  turbulence  with  good  accuracy  during  daytime. 

RESUME 

Nous  presentons  un  modele  simple  de  la  constante  de  structure  des  variations  de  l'indice  de 
refraction  atmospherique  (Cn2),  de  taille  caracteristique  ( lQ )  et  applicable  a  quelques 
centaines  de  metres  au-dessus  du  sol.  Les  donnees  d’entree  du  modele  sont  des  parametres 
meteorologiques  standard,  auxquels  on  ajoute  le  rayonnement  solaire  et  les  parametres 
caracteristiques  du  terrain.  On  compare  les  resultats  obtenus  aux  mesures  de  Cn2  et  de  l0 
effectuees  h.  l'aide  d'un  scintillometre.  Les  resultats  montrent  que  le  modele  permet  de 
prevoir  les  turbulences  atmospheriques  avec  une  bonne  precision  durant  le  jour. 
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EXECUTIVE  SUMMARY 

The  performance  of  electro-optical  systems  used  for  surveillance  and  target  detection  can 
be  degraded  due  to  the  effects  of  atmospheric  turbulence.  The  turbulent  atmosphere  induces 
slight  fluctuations  in  the  direction  of  light  rays  propagating  over  long  distances.  These  fluc¬ 
tuations  result  in  image  distortion  and  blur  and  may  limit  the  angular  resolution  of  an 
imaging  system.  The  effect  of  atmospheric  turbulence  is  often  seen  as  “shimmering”  of 
distant  objects  seen  through  binoculars  on  hot  sunny  days.  Atmospheric  turbulence  also 
affects  laser  beams,  inducing  intensity  scintillation  and  beam  wander,  effects  that  can 
degrade  the  performance  of  target  designators  and  beam  riders. 

The  intensity  of  atmospheric  turbulence  is  characterized  by  the  refractive  index  structure 
parameter  Cn  which  is  essentially  the  standard  deviation  of  differences  in  the  refractive 
index  of  air  measured  simultaneously  at  two  points.  This  quantity,  along  with  the  turbu¬ 
lence  inner  scale  lQ,  are  sufficient  to  calculate  the  optical  effects  of  military  significance. 
This  report  describes  a  model  that  predicts  Cn2  and  iQ  as  a  function  of  ordinary  meteorolog¬ 
ical  parameters,  solar  irradiance  and  some  properties  characteristic  of  the  local  terrain. 

The  impact  of  this  work  is  that  atmospheric  turbulence  can  be  predicted  with  a  simple 
model  for  situations  of  relevance  to  the  Canadian  Land  Forces.  Such  a  model  can  be  used 
in  simulations  and  studies  of  the  effectiveness  of  EO  weapons  and  surveillance  systems.  As 
the  angular  resolution  of  EO  systems  increases,  atmospheric  turbulence  is  more  frequently 
the  limiting  factor  in  system  performance.  Furthermore,  advances  in  adaptive  optics  tech¬ 
niques,  currently  used  in  astronomy  to  dynamically  correct  images  for  turbulence  effects, 
will  eventually  provide  technology  suitable  for  land  surveillance  applications.  This  report 
represents  a  First  step  in  understanding  atmospheric  turbulence  and  its  effects  on  EO 
systems. 
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1.0  INTRODUCTION 

Electro-optical  (EO)  systems  depend  on  the  propagation  of  light  over  large  distances  and 
are  thus  susceptible  to  the  influence  of  optical  turbulence.  The  turbulent  atmosphere  causes 
the  intensity  of  a  light  beam  to  fluctuate  or  scintillate,  causes  beam  wander  and  causes  the 
distortion  and  random  displacement  of  images.  Turbulence  can  limit  the  angular  resolution 
of  an  imaging  system  to  values  of  hundreds  of  microradians,  many  times  the  diffraction 
limit  of  visible  imagers,  thus  resulting  in  significantly  degraded  performance  (Ref.  1). 
Predictions  of  the  refractive  index  stmcture  parameter  Cn2  and  the  inner  scale  the  most 
important  parameters  which  characterize  the  optical  effects  of  atmospheric  turbulence,  can 
be  used  to  predict  and  optimize  the  performance  of  EO  systems. 

We  present  a  simple  model  of  atmospheric  turbulence  based  on  the  Monin-Obuhkov  simi¬ 
larity  theory  (Refs.  2  and  3)  as  interpreted  by  Thiermann  and  colleagues  in  Refs.  4  to  6.  The 
model  provides  height  profiles  of  Cn2  and  lQ  using  standard  meteorological  parameters  as 
inputs.  Model  results  are  compared  with  measurements  of  Cn2  and  lQ  made  with  a 
displaced-beam  scintillometer  installed  in  a  field.  The  scintillometer  provided  the  path 
averaged  values  of  Cn2  and  iQ  over  a  185  m  path  at  a  height  of  1.8  m.  Agreement  between 
the  model  and  measurements  are  generally  quite  good  under  typical  daytime  conditions.  At 
night,  under  stable  conditions  with  low  wind  speed,  the  model  exhibits  the  instability 
inherent  in  the  Monin-Obuhkov  similarity  theory. 

This  work  was  carried  out  at  DREV  from  Jan.  1994  to  Dec.  1995  under  PSC  32E15  “EO 
Battlefield  Surveillance/ Atmospheric  Propagation  Studies”. 

2.0  OPTICAL  TURBULENCE 

When  light  propagates  through  turbulent  air,  it  encounters  inhomogeneities  in  the  refractive 
index  that  transform  its  phase  front.  The  inhomogeneities  are  caused  by  the  mixing  of  air 
of  different  temperature  by  wind  and  convection.  Near  the  ground,  the  intensity  of  temper¬ 
ature  fluctuations  depends  on  the  type  of  ground  cover,  ambient  temperature,  solar  irradia¬ 
tion  and  wind  speed. 

Atmospheric  turbulence  may  be  characterized  by  three  parameters,  the  outer  and  inner 
scales  L0  and  £Q  and  the  refractive  index  structure  parameter  Cn2.  The  outer  scale  L0  repre¬ 
sents  the  minimum  scale  size  over  which  turbulence  energy  is  injected  into  the  air  through 
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convection  or  wind  shear.  Typical  values  for  L0  are  from  meters  to  hundreds  of  meters. 
However,  in  the  first  few  meters  above  the  ground,  the  surface  layer,  Lo  is  given  approxi¬ 
mately  by  L0=0.4h  where  h  is  the  height  above  the  ground  (Ref.  1).  The  continuous  motion 
of  turbulent  air  causes  the  turbulence  cells  to  continuously  break  up  into  ever  smaller  cells 
until  a  minimum  size  lQ  is  reached.  tQ  marks  the  turbulence  scale  size  below  which  viscous 
dissipation  converts  the  turbulent  energy  into  heat.  In  the  surface  layer  i0  is  of  the  order  of 
millimeters.  The  region  between  t0  and  L0  is  known  as  the  inertial  subrange. 


Cn2  is  a  measure  of  the  intensity  of  the  refractive  index  fluctuations.  For  homogeneous  and 
isotropic  turbulence  it  is  defined  as  (Ref.  1) 


C 


2 

n 


<  [n  (r j )  -n(r2)]2) 


(1) 


where  n  is  the  index  of  refraction  at  points  ?  j  and  r2 .  In  SI  units  (used  through  this  paper), 
Cn2  has  the  dimensions  of  m'2/3  and  varies  from  about  10" 16  to  10" 12  m”2/3  in  the  atmo¬ 
spheric  boundary  layer.  Equation  1  is  valid  for  distances  |?j  -  ?2|  much  greater  than  iQ  and 
much  smaller  than  L0. 

Cn2  is  expressed  in  terms  of  the  fluctuations  in  air  temperature  by  (Ref.  1) 
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where  CT2  (K2m'2/3)  is  the  temperature  structure  parameter,  P  is  the  air  pressure  in  hPa  (or 
millibars)  and  T  is  the  ambient  temperature  in  degrees  K. 


The  inner  scale  t0  (m)  is  give  by  (Ref.  6): 


,  _  .  3/4  -1/4 

l  =  7.4v  £ 

o 


(3) 


where  v  is  the  kinematic  viscosity  of  air  (1.46x  10"5  m2/s  at  STP)  and  £  (m2/s3)  is  the  dissi¬ 
pation  rate  of  turbulent  kinetic  energy.  CT2  and  £  are  provided  by  the  model  and  Eqs.  2 
and  3  are  then  used  to  calculate  Cn  and  6Q. 
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3.0  TURBULENCE  MODEL 

The  solutions  for  Cx2  and  £  are  expressed  in  terms  of  the  friction  velocity  u*  (m/s),  the 
turbulent  temperature  scale  T*  (K)  and  the  Monin-Obuhkov  length  L  (m).  u*,  T*  and  L  are 
scaling  parameters  defined  as  (Ref.  6): 


- 

(z  } 

- 

u*  =  uk 
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-V(L) 
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\Zo) 

_ 
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(4) 
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(6) 


In  Eqs.  4  to  6,  u  (m/s)  is  the  wind  velocity  measured  at  height  zu  (m)  above  the  ground,  zG 
(m)  is  the  roughness  length  of  the  ground  surface,  Q0  (Km/s)  is  the  vertical  turbulent  kine¬ 
matic  heat  flux,  T  is  the  air  temperature  measured  at  height  zu,  k  (dimensionless)  is  the  Von 
Karman  constant  taken  to  be  equal  to  0.35  and  g  is  the  acceleration  of  gravity  (9.81  m/s  ). 


Derivation  of  the  scaling  parameters  depends  upon  the  form  of  xjr(L)  which  is  chosen 
according  to  whether  the  conditions  are  stable  (air  warmer  than  the  ground)  or  unstable  (air 
colder  than  the  ground).  The  \y(L)  are  given  by  (Refs.  5  and  7): 


V(L)  =  -4.7zu/L 


y(L)  =  2  In 
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zu/L  >  0,  stable  (7) 
zu/L  <  0,  unstable  (8) 


where  y  = 


(l-15zu/L) 


1/4 


During  the  day  Q0  is  estimated  by  (Refs.  5  and  8) 


Q0  = 


JL 

CpP 
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(9) 
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where  rj  (dimensionless)  is  an  empirical  constant  equal  to  0.9,  cp  (J/kg/K)  is  the  specific 
heat  of  air,  p  (kg/m3)  is  the  density  of  air,  a  (dimensionless)  is  a  parameter  between  0  and 
1  that  represents  the  humidity  of  the  ground,  y  is  the  ratio  of  the  specific  heat  of  air  at 
constant  pressure  to  the  latent  heat  of  water  vapor,  s  is  the  temperature  derivative  of  satu¬ 
ration  specific  humidity.  The  dimensionless  ratio  y/s  is  a  function  of  temperature  only  and 
is  tabulated  in  Ref.  8.  A  (dimensionless)  is  the  surface  albedo,  R  (W/m2)  is  the  solar  irradi- 
ance  and  (3  is  an  empirical  constant  which  we  have  set  equal  to  35  W/m  . 

Equation  9  describes  how  the  solar  irradiation  is  partitioned  at  the  ground  leaving  a  compo¬ 
nent  that  contributes  to  the  sensible  heat  flux.  The  factor  (l-A)R  is  the  fraction  of  solar  irra¬ 
diation  absorbed  by  the  ground.  Part  of  this  power  is  dissipated  by  infrared  radiation  and 
ground  conduction.  This  loss  is  accounted  for  by  the  factor  r|  and  offset  p.  The  remaining 
power  gives  rise  to  both  sensible  and  latent  heat  fluxes.  The  latent  heat  flux  is  due  to  both 
evaporation  at  the  ground  surface  and  transpiration  by  vegetation.  The  fraction  of  power 
available  for  sensible  heat  flux  is  determined  by  the  factor  l-oc/(l+y/s),  plotted  in  Fig.  1. 
According  to  Holtslag  and  Van  Ulden  (Ref.  8),  a~l  is  an  appropriate  value  for  grass 
covered  land  as  is  the  case  at  our  location.  For  dry  areas  without  vegetation  (desert),  a~0. 
For  a=l  at  a  temperature  of  25°  C,  l-a/(l+y/s)=0.26  which  means  that  only  just  over  one 
quarter  of  the  available  solar  irradiance  is  transformed  into  sensible  heat  flux. 


FIGURE  1  Fraction  of  power  available  for  sensible  heat  flux  as  a  result  of  evaporation 
loss,  a  is  the  ground  humidity  parameter. 
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The  nighttime  kinematic  heat  flux  is  difficult  to  estimate  because  it  does  not  depend  on  one 
dominant  factor  as  is  the  case  with  solar  irradiance  for  the  daytime  heat  flux.  The 
phenomena  on  which  the  nighttime  heat  flux  depends;  radiative  cooling  of  the  ground, 
thermal  conduction  losses  and  latent  cooling  or  heating  due  to  evaporation  or  formation  of 
dew,  may  all  have  similar  importance.  Not  only  is  the  interaction  of  these  parameters 
complex  but  information  regarding  them  would  be  too  specialized  for  a  practical  model.  To 
estimate  the  nighttime  heat  flux,  Thiermann  and  Kohnle  (Ref.  5)  developed  an  expression 
which  depends  on  the  one  parameter  that  is  both  important  and  readily  available;  the  wind 
speed  u.  The  expression  is 


1  + 


3 

cu 


c  3 

Q  U 

xomax 


(10) 


where 


c  = 


k2T 


27  5gzu(ln(zu/z0))2 


The  structure  of  Eq.  10  ensures  that  Q0  does  not  diverge  at  high  wind  speed  but  increases 
asymptotically  to  a  saturation  value  Qomax.  This  saturation  value  is  given  by  Qomax 
=Hmax7cpP  where  the  sensible  heat  flux  Hmax  is  chosen  to  be  between  -5  to  -100  W/m  .  The 
most  appropriate  value  of  Hmax  depends  on  local  vegetation,  ground  characteristics,  air 
humidity  and  cloud  cover.  For  the  calculations  in  this  paper  we  set  Hmax  =-12  W/m  ,  but 
as  will  be  shown  later,  Q0  given  by  Eq.  10  was  in  poor  agreement  with  Q0  measured  with 
the  scintillometer  regardless  of  the  value  of  Hmax. 

Since  the  Monin-Obuhkov  length  L  cannot  be  written  as  an  explicit  function  of  \| /,  an  iter¬ 
ative  procedure  must  be  used  to  find  mutually  consistent  solutions,  u*  and  L,  to  Eqs.  4  to  6. 
Following  Holstag  and  Van  Ulden  (Ref.  8),  an  initial  guess  for  u*  is  found  by  evaluating 
Eq.  4  with  vjr(L)=0.  L  is  then  estimated  by  using  the  initial  u*  in  Eq.  6.  u*  is  then  reevaluated 
using  L  with  \|/(L)  given  by  Eqs.  7  or  8,  thus  beginning  a  cycle  where  subsequent  values  of 
u*  are  compared  to  the  previous  one  until  they  are  within  0.5%  of  each  other.  Typically 
three  to  ten  iterations  are  required  to  determine  u*  and  L. 
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The  decision  as  to  whether  conditions  are  stable  or  unstable  is  based  on  the  value  of  the 
daytime  heat  flux  Q0  given  by  Eq.  9.  If  the  daytime  Q0  is  positive  then  conditions  are 
assumed  to  be  unstable  and  Eq.  8  is  used  to  derive  L.  If  the  daytime  Q0  is  negative  then 
conditions  are  assumed  to  be  stable,  the  night  time  heat  flux  is  calculated  using  Eq.  10,  and 
Eq.  7  is  used  to  calculate  L. 


Once  the  scaling  parameters  are  determined,  the  vertical  profiles  of  CT2  and  e  are  given  by 
(Ref.  6): 


CZ.  =  4p- 
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for  the  stable  case  (z/L>0)  and  by 
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(14) 


for  the  unstable  case  (z/L<0).  the  parameters  Cn2  and  lQ  of  interest  here  are  then  found 
substituting  CT2  and  £  into  Eqs.  2  and  3. 

A  multitude  of  physical  quantities  are  required  in  the  calculations.  The  specific  heat  of  air 
cp  is  taken  to  be  a  constant  equal  to  1004.6  J/kg/K  (Ref.  9).  The  density  of  dry  air  p,  the 
ratio  y/s  and  the  kinematic  viscosity  of  air  v  are  given  in  following  formulae  where  temper¬ 
ature  T  is  in  degrees  C: 
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p  (kg/m  )  =  1.286- 0.00405 T 


y/s  =  1.4631  -0.0923T  +  0.0027T2-3. 18x10  5T3 


% , ,  2.  ,  f  1.718  +  0.0049T^  1A-5 
v(m  /s)  =  -  xlO 


(15) 

(16) 

(17) 


Equation  15  is  derived  from  Annex  C  of  Ref.  9,  Eq.  16  is  a  least  squares  fit  to  data  tabulated 
in  Ref.  8  and  Eq.  17  is  from  Annex  C  of  Ref.  9.  The  roughness  length  zQ  chosen  to  represent 
our  location  was  0.02  m  based  on  estimated  values  listed  in  Ref.  8.  The  ground  humidity 
parameter  a  was  set  at  1 .  The  value  of  ground  albedo  A  used  is  0.2  based  on  measurements 
made  with  a  commercial  albedometer  and  is  a  typical  value  for  grass  covered  land 
according  to  Ref.  8.  For  the  purposes  of  comparing  the  model  to  measured  turbulence,  we 
have  used  the  measured  solar  irradiance  R  as  an  input  to  the  model.  Flowever  R  may  be 
easily  estimated  on  the  basis  of  geographic  location,  time  of  day  and  cloud  cover.  An  algo¬ 
rithm  to  calculate  R  is  given  in  Annex  A. 


4.0  TURBULENCE  MEASUREMENTS 

The  experiments  were  performed  at  the  Defence  Research  Establishment  Valcartier 
(DREV)  in  July  and  August,  1995.  The  site,  known  as  Parc  Lemay,  is  a  level  grass  covered 
field  with  sandy  soil.  It  is  about  300  m  across,  surrounded  by  spruce  trees  on  the  West  and 
North  sides  and  by  low  buildings  on  the  East.  The  scintillometer  transmitter  and  receiver 
were  housed  in  small  instrument  sheds  with  a  propagation  path  of  1 85  m.  The  path  was 
1.8  m  above  the  ground  and  oriented  approximately  North-South  with  prevailing  winds 
coming  from  the  West.  The  layout  of  the  experiment  is  shown  in  Fig.  2. 

The  Scintec  SLS-20  scintillometer,  used  to  measure  Cn2  and  6Q,  is  described  in  detail  in 
Refs.  10  and  11.  The  SLS-20  uses  the  independently  measured  scintillation  of  two 
displaced  but  parallel  laser  beams  to  derive  the  turbulence  parameters.  The  source  is  a 
1  mW  laser  diode  with  a  wavelength  of  670  nm  and  divergence  of  approximately  5  mrad. 
The  beam  displacement  of  2.7  mm  is  achieved  by  passing  the  beam  through  a  birefringent 
calcite  crystal  which  splits  the  beam  into  displaced  but  parallel  components  with  crossed 
polarizations.  At  the  receiver  the  two  overlapping  beams  are  resolved  by  a  polarizing  beam 
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splitter,  allowing  the  scintillation  of  each  beam  to  be  measured  independently  by  two 
2.5  mm  diameter  detectors.  The  transmitted  laser  beams  are  modulated  at  20  kHz  to  allow 
AC  coupling  of  the  detector  outputs  which  eliminates  the  background.  The  calculations  to 
derive  Cn2  and  t0  are  carried  out  in  real  time  by  the  instruments  laptop  computer  which  also 
derives  CT2,  the  kinetic  energy  dissipation  rate,  the  turbulent  fluxes  of  sensible  heat  and 
momentum,  and  the  Monin-Obukhov  length  using  the  expressions  of  Sec.  2.  A  calibration 
was  carried  out  for  the  1 85  m  path  length  according  to  the  manufacturer’s  instructions.  The 
scintillometer  receiver  is  shown  in  Fig.  3. 


FIGURE  2  Plan  of  Parc  Lemay,  DREV,  where  turbulence  measurements  were  made. 
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FIGURE  3  Receiver  of  Scintec  SLS-20  scintillometer  receiver  installed  for 
measurements. 


FIGURE  4  The  anemometer  and  the  pyrradiometer  (for  measuring  surface  albedo) 
installed  near  the  midpoint  of  the  scintillometer  path. 
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Meteorological  conditions  were  monitored  simultaneously  near  the  scintillometer  path 
using  standard  sensors  for  wind  speed,  wind  direction,  barometric  pressure,  solar  irradia¬ 
tion,  air  temperature  and  dew  point  temperature.  The  ground  albedo  was  also  measured 
with  a  pyrradiometer  or  “albedometer”.  Some  of  the  instruments  and  the  character  of  the 
experiment  site  are  shown  in  Fig.  4.  The  wind  measurements  were  made  at  a  height  of  3.0 
m.  The  temperature  and  dew  point  measurements  were  made  at  a  height  of  1 .5  m.  The  mete¬ 
orological  data  was  recorded  using  the  SLS-20  analog  input  channels.  The  meteorological 
sensors  used  are  summarized  in  Table  I. 


Table  I 

Meteorological  sensors  used  in  experiment 


Physical  quantity 

Instrument  name 

Manufacturer 

Model 

Accuracy 

A  time 
(s) 

Wind  speed  & 
direction 

Wind  monitor 

Young 

05103 

±0.25% 

not 

specified 

Temperature  & 
dewpoint 

Dewpoint  & 
temp,  monitor 

General  Eastern 

1200  MPS 

±  0.2°  C 

1.6  °C/s 

Atm.  pressure 

Barometer 

Qualimetrics 

7105-A 

±  1% 

not 

specified 

Rain  rate 

Opt.  precip.  sens. 

STI 

ORG-705 

±  4% 

10 

Pyranometer 

Matrix 

R-1950 

±  5% 

0.001 

Ground  albedo 

Pyrradiometer 

Qualimetrics 

3040-A 

not 

specified 

40 

5.0  COMPARISON  OF  MODEL  AND  MEASUREMENTS 

Figures  6  shows  measurements  during  a  four  day  period  from  midnight  local  time  Aug.  23 
to  midnight  Aug.  27,  1995.  The  first  day  was  overcast  with  subsequent  days  mostly  clear 
as  indicated  by  the  solar  irradiance  curve  (a).  Wind  speed  (b)  and  relative  humidity  (c) 
exhibit  diurnal  cycles  with  wind  usually  less  than  3  m/s  at  night  and  between  3  and  8  m/s 
during  the  day.  Humidity  was  often  100%  at  night  and  always  less  than  75%  during  the  day. 
No  record  was  kept  of  the  cloud  conditions  at  night. 
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FIGURE  5  Measurements  from  midnight  local  time  Aug.  23  to  midnight  Aug.  27,  1995: 

solar  irradiance  (a),  wind  speed  at  3  m  height  (b)  and  one  hour  averages  of  rela¬ 
tive  humidity  (c). 

Measurements  of  Cn2  and  model  calculations  are  shown  in  Figs.  6  (a)  and  (b).  The  calculated 
Cn2  follows  closely  the  measured  data  for  daytime  periods.  At  night,  the  agreement  between 
measured  and  modelled  Cn2,  although  reasonably  good,  is  poorer  than  during  the  day  due 
primarily  to  the  approximate  method  of  determining  the  nighttime  heat  flux  based  on  the  wind 
speed  (Eq.  10).  Gaps  in  the  nighttime  model  results  occur  when  the  wind  speed  is  less  than 
0.5  m/s.  Below  this  value,  calculations  were  not  performed  to  avoid  inherent  instabilities  in  the 
Monin-Obukhov  similarity  theory. 


cn2  (m-2'3)  Cn2  (m^3) 
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Measurements  and  model  calculations  of  6Q  are  shown  in  Figs.  7  (a)  and  (b).  The  relatively 
large  values  of  60  during  calm  night  periods  and  small  values  during  windy  day  conditions 
are  typical.  Agreement  between  the  measurements  and  model  is  generally  quite  good,  except 
at  night  when  the  wind  speed  is  less  than  0.5  m/s  and  calculations  were  not  performed.  That 
the  model  is  better  able  to  represent  nighttime  i0  than  Cn2  is  consistent  with  an  analysis  by 
Thiermann  and  Kohnle  (Ref.  5)  of  the  sensitivity  of  the  similarity  equations  to  input  errors. 
They  showed  that  Cj2,  and  hence  Cn2,  is  ten  times  more  sensitive  to  errors  in  input  Q0  than 

4  is- 

0.05 


0.04 

E 

0.03 

o 

O 

CD 

■|  0.02 

JZ 

05 

‘c 

0.01 


0 

01234567801  2345678 

wind  speed  (m/s)  wind  speed  (m/s) 

FIGURE  8  Measured  nocturnal  kinematic  heat  flux  as  a  function  of  wind  speed. 

Equation  10  is  plotted  for  several  values  of  saturation  heat  flux  Homax. 

The  main  reason  that  both  the  modelled  turbulence  parameters  are  in  poorer  agreement  with 
the  measurements  at  night  is  because  the  nighttime  heat  flux  estimated  with  Eq.  10  is  not  very 
accurate.  Figure  8  (a)  shows  the  kinematic  heat  flux  measured  with  the  scintillometer  as  a 
function  of  wind  speed.  The  curves  are  Eq.  10  plotted  for  various  values  of  saturation  heat 
flux  Homax.  For  our  calculations  we  set  Homax=-12  W/m2  although  it  is  clear  from  the  figure 
that  agreement  between  the  measurements  and  Eq.  10  is  poor  regardless  of  the  choice  of 
Homax-  The  scatter  of  the  data  points  shows  that  the  nocturnal  Q0  has  only  a  vague  wind  speed 
dependence  under  the  conditions  of  our  measurements. 
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The  "measured"  Q0  shown  in  Fig.  8  are  provided  by  the  Scintec  scintillometer  software  and 
are  derived  from  the  measured  Cn2  and  l0  via  Eqs.  11  to  14.  Therefore,  if  we  choose  the 
same  value  of  Q0  as  determined  by  the  scintillometer,  the  model  should  produce  the  same 
values  of  Cn2  and  lQ  as  recorded  by  the  scintillometer. 

It  is  possible  that  humidity  effects  play  a  role  in  generating  heat  flux,  something  which  is 
not  taken  into  account  in  Eq.  10.  The  measurements  indicate  that  periods  of  high  Cn  at 
night  coincided  with  100%  relative  humidity,  for  example  hours  48  to  54.  Dew  was  usually 
observed  in  the  mornings,  evidence  that  latent  heat  was  released  during  the  night. 
Figure  8(b)  shows  the  nocturnal  Q0  as  a  function  of  wind  speed,  but  with  data  removed  that 
corresponded  with  relative  humidity  greater  than  90%.  If  formation  of  dew  contributed  to 
the  measured  Q0  then  we  would  expect  that  the  wind  speed  dependence  be  more 
pronounced  for  data  restricted  to  relatively  dry  air.  However,  the  Q0  data  of  Fig.  8  (b)  is  not 
significantly  better  correlated  with  wind  speed  than  the  data  of  Fig.  8  (a). 

Further  work  is  thus  required  to  determine  the  role  (if  any)  of  humidity  changes  and  other 
possible  effects,  in  generating  nocturnal  heat  flux.  An  instrument  to  measure  Q0  indepen¬ 
dently  of  the  scintillometer  would  be  useful  in  this  regard.  Estimating  the  night  time  heat 
flux  is  the  weakest  part  of  the  model. 


6.0  CONCLUSIONS  AND  FUTURE  WORK 

We  presented  results  from  a  simple  model  for  the  structure  constant  of  atmospheric  refrac¬ 
tive  index  fluctuations  (Cn2)  and  inner  scale  (i0)  valid  in  the  lowest  few  hundred  meters 
above  land.  The  model  is  based  on  the  Monin-Obuhkov  similarity  theory  (Refs.  2  and  3). 
The  model  equations  follow  Thiermann  et.  al.  (Refs.  4  to  6).  The  required  input  parameters 
are  solar  irradiance,  air  temperature,  wind  speed,  barometric  pressure  and  terrain  character¬ 
istics;  albedo,  roughness  length  and  ground  humidity. 

Values  of  Cn2  and  iQ  calculated  with  the  model  were  in  good  agreement  with  measurements 
made  with  a  displaced-beam  scintillometer  installed  in  a  field.  The  scintillometer  provided 
the  path  averaged  values  of  Cn2  and  l0  over  a  185  m  path  at  a  height  of  1.8  m.  The  poorest 
agreement  was  observed  during  calm  night  periods  where  the  model  failed  to  predict  the 
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high  levels  of  intermittent  Cn2  measured  by  the  scintillometer.  It  is  believed  that  the  release 
of  latent  heat  at  night  during  formation  of  dew  may  contribute  to  the  nocturnal  heat  flux,  an 
effect  that  was  not  taken  into  account.  Also,  the  Monin-Obuhkov  similarity  theory  is  inher¬ 
ently  unstable  for  low  wind  speed  and  stable  night  time  conditions. 

Predictions  of  Cn  and  60  made  with  this  model  can  be  used  to  predict  and  optimize  the 
performance  of  electro-optical  systems  which  are  susceptible  to  the  scintillation,  beam 
wander  and  image  distortion  caused  by  optical  turbulence  over  land.  Our  next  effort  will  be 
to  include  humidity  fluctuations  in  the  model  which  will  allow  Cn2  to  be  calculated  in  the 
infrared  and  microwave  bands.  Also,  path-averaged  effects  on  EO  systems  such  as  scintil¬ 
lation,  point-spread  function  and  image  displacement  will  be  included  explicitly  in  the 
model  and  compared  to  measurements  of  image  quality  presently  being  made  at  the  Parc 
Lemay  test  site. 
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ANNEX  A 

CALCULATION  OF  GROUND  LEVEL  SOLAR  IRRADIANCE 

The  solar  irradiation  R  at  the  ground  can  be  estimated  for  a  given  time,  date  and  geograph¬ 
ical  location  by  calculating  the  solar  elevation  angle  9  which  is  the  angle  of  the  sun  above 
the  local  horizon.  The  downwelling  radiation  at  the  surface  is  then  given  by: 


R  =  ST  sin  0  for  daytime  (0  positive)  (Al) 

R  =  0  for  nighttime  (9  negative) 


where  S  is  the  solar  constant  or  the  intensity  of  solar  radiation  at  the  top  of  the  atmosphere 
and  T  is  the  sky  transmissivity.  S  will  be  taken  to  be  1370  W/m  and  T  is  given  approxi¬ 
mately  by  (Ref.  9): 

T  =  0.6 +  0.2  sinG  .  (A2) 


From  Ref.  9,  the  expression  for  the  local  solar  elevation  angle  is 


sinG  =  sin (J) sin 5  -  cos (j) cos 8  cos 


TCt 


utc 


12 


(A3) 


where  (j)  and  a  are  the  latitude  (positive  north)  and  longitude  (positive  west)  in  radians.  5  is 
the  solar  declination  angle  or  the  angle  of  the  sun  above  the  equator,  in  radians  and  tutc  is 
the  Coordinated  Universal  Time  in  hours.  5  is  given  by: 


8 


<|>rCOS 


f2TC(d-dr) 


(A4) 


where  (j)r  is  the  latitude  of  the  Tropic  of  Cancer  (23.45°  or  0.409  radians)  d  is  the  number 
of  the  day  of  the  year  (e.  g.  August  27  is  day  number  239)  and  dr  is  the  day  of  the  summer 
solstice,  173.  dy  is  the  average  number  of  days  in  a  year,  365.25. 


Combining  the  above  expressions  and  using  the  geographical  location  of  DREV, 
<()=46053’5”  (0.818  rad)  and  a  =71°28’47”  (1.248  rad),  yields  the  solar  irradiance  curve 
for  August  27  shown  in  Fig.  Al.  The  data  is  plotted  as  a  function  of  local  Eastern  Daylight 
Savings  time  by  adding  3.75  hours  to  tutc. 
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FIGURE  A1  Example  of  calculated  and  measured  solar  irradiance  at  DREV. 

August  27,  1995  was  a  mostly  clear  day  with  some  thin,  high  cloud  in  the  morning.  It  can 
be  seen  from  the  figure  that  the  calculated  solar  irradiance  is  in  good  agreement  with  the 
measured  values. 

The  simple  model  presented  here  does  not  take  into  account  the  effect  of  cloud  cover.  A 
more  elaborate  expression  for  atmospheric  transmittance  T  (Eq.  A2)  exists  but  the  height 
and  thickness  of  the  clouds  must  be  estimated  which  is  too  unreliable  for  our  purposes. 
More  detail  on  this  subject  is  given  in  Ref.  9,  p.  257.  Also,  the  model  does  not  explicitly 
take  into  account  the  spectral  dependence  of  atmospheric  absorption.  The  semi-empirical 
expression  for  T,  Eq.  A2,  applies  to  the  “short  wave”  solar  radiation  in  the  spectral  band 
approximately  0.4  to  2  |im. 
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